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Abstract 

We study the problem of the Landau gauge fixing in the case of the SU{2) 
lattice gauge theory. We show that the probability to find a lattice Gribov copy 
increases considerably when the physical size of the lattice exceeds some critical value 
« 2.75/v^, almost independent of the lattice spacing. The impact of the choice of 
the copy on Green functions is presented. We confirm that the ghost propagator 
depends on the choice of the copy, this dependence decreasing for increasing volumes 
above the critical one. The gluon propagator as well as the gluonic three-point 
functions are insensitive to choice of the copy (within present statistical errors). 
Finally we show that gauge copies which have the same value of the minimisation 
functional (/ d'^x{A'^)'^) are equivalent, up to a global gauge transformation, and 
yield the same Green functions. 
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1 Introduction 

The question of the quantisation in the infrared of a non-Abelian gauge theory is very 
important for the understanding of diverse non-perturbative phenomena. The main prob- 
lem that arises v^^hile performing the quantisation in the covariant gauge is the problem of 
the non-uniqueness of the solution of the equation specifying the gauge-fixing condition, 
that was pointed out by Gribov |T|. In fact, the solution is unique for fields having small 
magnitude. A solution proposed in [l] consists in restricting the (functional ) integration 
over gauge fields in the partition function to the region inside the so called Gribov horizon 
(the surface, closest to the trivial field A^j, = 0, on vv^hich the Faddeev-Popov determinant 
vanishes). 

The Gribov's quantisation prescription is explicitly realised in the lattice formulation. 
Indeed, the procedure of Landau gauge fixing on the lattice (see [2] for a revievv^) con- 
straints all the eigenvalues of the Faddeev-Popov operator to be positive. Thus lattice 
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gauge configurations in Landau gauge are located inside the Gribov region, and all calcu- 
lations on the lattice are done within the Gribov's quantisation prescription. But we also 
know that this is not enough - the gauge is not uniquely fixed even in this case, and there 
are secondary solutions (Gribov copies) inside the Gribov region The smaller 

region free of Gribov copies is called the fundamental modular region (0,0, On a 
finite lattice the gauge may be fixed in a unique way, but that reveals to be technically 
difficult (cf. the following section). 

The quantities that are sensitive to the choice of the infrared quantisation prescription 
are Green functions. A lot of work regarding the infiuence of Gribov copies on lattice ghost 
and gluon propagators in the case of SU (2) and SU (3) gauge groups has already been done 
(L7j,[8j,[9j,[10j,[llJ,[r2],[ini)- Summarising the results, the infrared divergence of the ghost 
propagator is lessened when choosing Gribov copies closer to the fundamental modular 
region, and the gluon propagator remains the same (within today's statistical precision). 
In this paper we confirm these results in the case of the SU (2) gauge group, and present 
a study of the infiuence of Gribov copies on the three-gluon vertex in symmetric and 
asymmetric kinematic configurations. We also discuss in details the structure of minima 
of the gauge-fixing functional. The authors of [H] reported that the number of minima of 
the gauge-fixing functional decrease with the [3 parameter (the bare lattice coupling). We 
perform a thorough analysis of the volume dependence of this phenomenon. We define a 
probability to find a lattice Gribov copy, and our computation shows that this probability 
increases considerably when the physical size of the lattice exceeds some critical value, 
around 2.75/v^, where a is the string tension. 



2 Non-perturbative Landau gauge fixing on the lattice 

A lattice gauge configuration Ucq generated during the simulation process is not gauge- 
fixed. One has to perform a gauge transformation {m(x)} on it in order to move it along its 
gauge orbit up to the intersection with the surface fi \U^{x)] = specifying the Landau 
gauge-fixing condition. But there is no need to have an explicit form of {u{x)}. Instead 
we do an iterative minimisation process that starts at Ucq and converges to the gauge 
fixed configuration U^^^ Let us first illustrate it on the example of the Landau gauge in 
the continuum limit. For every gauge field A one defines a functional 

F4[m(x)] = -Tr J d''xA^;\x)A^;\x) = -||^("^||', uix) e SUiN,). (1) 

Expanding it up to the second order around some group elements {^0(2;)} we have (writing 
u = e^uo, X G su{Nc)) 



(2) 



Note that the quadratic form defining the second order derivatives is the Faddeev-Popov 
operator. Obviously, if Uq is a local minimum we have a double condition 

/l("0) _ Q 

A^FP [^^"°^] is positive definite. ^ ' 

Hence, the minimisation of the functional ^ allows not only to fix the Landau gauge, 
but also to obtain a gauge configuration inside the Gribov horizon (cf. Figure HI). On the 
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Figure 1 : Restriction of the integration domain in the partition function to the Gribov region (hatched) , 
and the Gribov horizon. "FP" is the Faddeev-Popov operator. 



lattice, the discretised functional fH) reads 

FuM = -^'y^tTrJ2u{x)U^{x)u\x + e^), (4) 

X,fJ, 

where f/^(x) is the standard gauge link variable. Then at a local minimum uq we have a 
discretised Landau gauge fixing condition 

that we write in a compact form V^^^T"'' = 0. Here 



The second derivative (the equivalent of the second-order term in ^) can be written 



as 

^,FuHs,x)] = -^{u,V,A'^) (7) 



where 



A', = l [~u{x)U^;\x) + U^-\x)u:{x + e,) + u{x + e,)U^;^\x) - f7j"»(^)^(^)) • (8) 

This defines a quadratic form {oj , M^pp\U]u}) where the operator A^^p[[/] is the lattice 
version of the Faddeev-Popov operator It reads 

{m%[U]u)\x) = i {u'ix + e,) - u'ix)) - {x ^ x - e,)+ 
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where 



S;\x) = --Tv{{t^,t'}{U,ix) + Ulix))) 



(10) 



It is straightforward to check that in the continuum limit a — > one finds the familiar 
expression of the Faddeev-Popov operator. 

For the Landau gauge fixing in our numerical simulation we have used the Overrelax- 
ation algorithm with uj = 1.72. We stop the iteration process of the minimising algorithm 
when the following triple condition is fullfield: 



ViN^ - 1) 



5^Tr[(V,4-)) (V,4-))^ 
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where u^"\x) is the matrix of the gauge transformation u{x) at the iteration step n, and 
the charge 



A^(x,t) = / d^xA'^{x,t) 



(12) 



must be independent of t in Landau gauge when periodical boundary conditions for the 
gauge field are used. The choice of numerical values for the stopping parameters is dis- 
cussed at the end of the next section. 



3 The landscape of minima of the functional Fjj 

The functional Fu ^ has a form similar to the energy of a spin glass. The last is 
known to posses a very large number of metastable states, i.e. spin configurations whose 
energy increases when any spin is reversed. Typically, the number of these states grows 
exponentially with the number of spins. 

Let us consider the landscape of the functional Fjj. One of its characteristics is the 
distribution of values at minima F-^m of Fu. We know that for small magnitudes of the 
gauge field all the link matrices U^{x) (they play the role of couplings between the "spin" 
variables) are close to the identity matrix, and thus the minimum is unique. Their number 
increases when the bare lattice coupling /3 decreases, because the typical magnitude of the 
phase of U^{x) grows in this case and thus link matrices move farther from the identity 
matrix. The number of minima also increases with the number of links (at fixed /?) because 
in this case there are more degrees of freedom in the system. At FigureElwe present typical 
histograms of the distribution of Fmin for one given gauge configuration, as a function of 
the l3 parameter. We see that when /? is large (/5 > 2.5 for the volume = 8*^ considered 
at Figure 121), we typically find only one value of -Fmin- On the contrary, for very small 
values of (3 we find many different minima (cf. Figure [2]). These values of [3 correspond 
to almost random links U^{x). In the following section we show that field configurations 
having the same F-^m are in fact equivalent up to a global gauge transformation. 

The probability to find a secondary minimum depends on the value of the /3 parameter. 
We can calculate this probability in the following way. For each orbit we fix the gauge 
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Fmin distribution, SU(2) 8"^ 100 entries/hist 
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Figure 2: Histogram of minima values F^un for different /3 for one particular gauge-field configuration, 
|(9^A°(a;)| < 10~^. Here TVqf = 100. The histograms are ranged in the order the increase of the P 
parameter. 



Nqp times, each gauge fixing starts after a random gauge transformation of the initial 
field configuration. We thus obtain a distribution of minima -Fmin- This distribution gives 
us the number of minima A^(F*) as a function of the value of = F^^^. The relative 
frequency of a minimum is defined by 

where — ^gf- Then the weighted mean number of copies per value of Fmin is 

given by 

N = ^uJiN{F'). (14) 

i 

This allows us to define a probability to find a secondary minimum when fixing the Landau 
gauge for a gauge field configuration : 

Piconf = 1 - Y~N{F^y ^^^^ 

If one finds the same value of Fmin for all N^p tries then this probability is zero. On the 
contrary, if all Fj are different then p is close to one. 

Having the probability to find a secondary minimum when fixing the gauge for one 
particular configuration we can calculate the Monte-Carlo average {-k) on gauge orbits, i.e. 
on "spin couplings" U^{x). We obtain finally the overall probability to find a secondary 
minimum during a numerical simulation: 

W)-(Pico„f) = l-(^-|^). (16) 
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F^i„ distribution, SU(2) V=8* p=10"^ 250 entries 




Figure 3: Histogram of minima values i^min for very a small value of /3 = 10 < 10 



We have performed simulations in the case of SU{2) latice gauge theory at volumes 
V = {8^, 10"^, 12^, W^} for P varying from 1.4 to 2.9. For each value of /3 we generated 100 
independent Monte-Carlo gauge configurations, and we fixed the gauge Nqp = 100 times 
for every configuration. Between each gauge fixing a random gauge transformation of the 
initial gauge configuration was performed, and the minimising algorithm stops when the 
triple condition (fTT|l is satisfied. The resulting probability to find a secondary minimum 
is presented at (Figures SHU). 

As expected, the probability is small when P is large, and it is close to one when 
P is small. The dispersion was calculated using the Jackknife method. The physical 
meaning of this dispersion is the following: when the error is small, all gauge configurations 
have a similar number of secondary minima. On the contrary, this dispersion is large if 
there are some exceptional gauge configurations having a different number of copies. At 
small P almost all gauge configurations have many secondary minima, that is why the 
dispersion of the probability is small. At large P almost all gauge configurations have 
a unique minimum, but some of them can have copies. This may considerably increase 
the dispersion of the probability. The appearance of exceptional gauge configuration 
possessing a large density of close-to-zero eigenvalues of the Faddeev-Popov operator has 
been recently reported j^. Probably these fields are related with those having a lot of 
secondary minima at large /?'s, and this correlation deserves a separate study. 

We can fit the data from Figures SHU with an empirical formula 

^C) = TT^^ (17) 

in order to define a characteristic coupling Pc when the probability to find a copy decreases 
considerably. One can define Pc as corresponding to the semi-heights of the probability 
function P{P). At this value of P an equally probable secondary attractor of the functional 
Fu appears. The fit has been performed for the points between dashed lines at Figures 
mm and the results for the fit parameters are given ibidem and in Table [H We see that 
Pc depends on the volume of the lattice. Let us check whether these values correspond to 
some physical scale. According to works one has the following expression for the 
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Figure 4: Probability (averaged over gauge orbits) to find a secondary minimum as a function of f3 at 
volume V — 8^. 




Figure 5: Probability (averaged over gauge orbits) to find a secondary minimum as a function of /3 at 
volume V — 10'*. 
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SU(2) V=12^p,A^(x)|<10- 
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Figure 6: Probability (averaged over gauge orbits) to find a secondary minimum as a function of f3 at 
volume V — 12^. 



SU(2) V=16''|3,A^(x)|<10- 
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Figure 7: Probability (averaged over gauge orbits) to find a secondary minimum as a function of /3 at 
volume V — 16'*. 
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string tension cr for /? > 2.3: 



era 



with c = 4.38(9) and d = 1.66(4). Using this formula, we define a characteristic scale 
corresponding to the critical values /3c from (Figures IHEI) : 

Ac = a{Pc) ■ L 

in the string tension units, L is the length of the lattice. In the last column of Table [T] 
we summarise the results. We see that for the values of Pc in the scaling regime when the 



L 




xVndf 


Ac, in units of \l\fo 


8 


2.221(14) 


0.27 


3.85(31) 


10 


2.342(6) 


0.28 


3.20(21) 


12 


2.44(1) 


0.22 


2.78(20) 


16 


2.541(5) 


1.92 


2.68(17) 



Table 1: The characteristic length defining the appearance of secondary minima. The errors for Ac 
include errors for d and c parameters, and the fitted error for /3c 



formula (fT8|) is applicable (/3 > 2.3) we obtain compatible values for the physical length 
Ac. 

This suggests that lattice Gribov copies appear when the physical size of the lattice 
exceeds a critical value of around 2.75/i/o". At fist approximation Ac is scale invariant, 
but a slight dependence in the lattice spacing remains. 

In principle, the parameter /?c can be calculated with good precision. One should do 
it in the case of S'?7(3) gauge group, because the dependence of the lattice spacing on the 
bare coupling is softer, and the scaling of the theory has been better studied than in the 
case of the SU (2) theory. 

Let us discuss the dependence of above results on the choice of the stopping parameters 
(fTT| . In the Table El we present the probabilities P(2.0), P(2.3), P(2.8), calculated on a 
8^ lattice with different values of the parameters (fTI| . One can see that our choice of 
stopping parameters in (fTI|l is strict enough to ensure the independence of our results on 
the further increase of the parameter ©max^^ la^A^i- 



4 Green functions and the lattice Gribov copies 

4.1 Two- and three-point lattice Green functions 

Lattice Green functions are computed as Monte-Carlo averages over gauge-fixed gluon 
field configurations. 

Gluonic Green functions are defined using the definition Q of the gauge field. The 
ghost propagator F^'^^°'^[x — y) is computed, with the algorithm given in fTTj, by numerical 
inversion of the Faddeev-Popov operator Q, e.g. 

F^'^{x-yW'^l^M'fi\A\)'^^. (19) 
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The gluon propagator in Landau gauge may be parametrised as 

G^;}^\p, -P) ^ {Al{-p)Al{p)) = 5"^ (5,. - ^) (20) 
completed with 

G'g)'^^(0,0) = 5%,G(^)(0). (21) 

The ghost propagator is parametrised in the usual way: 

F'^'^^\p, -p) ^ (7^{-p)t{p)) = S^'F^'\p') (22) 

The coupling constant can be defined non-perturbatively by the amputation of a three- 
point Green-functions from its external propagators. But this requires to fix the kinematic 
configuration of the three-point Green- function at the normalisation point. On the lattice 
one usually uses either a fully symmetric kinematic configuration (denoted MOM) or 
a zero point kinematic configuration with one vanishing external momentum (denoted 

generically MOM). In what follows we only speak about gluonic three-point functions. 



Symmetric case 

There are only two independent tensors in Landau gauge in the case of the symmetric 
three-gluon Green function jH 



^j!U,M3(Pl'P2,P3) = 5miM2(Pi - P2)/., + 5m2M3(P2 - P3)mi + ^M3m(P3 " Pi) ^^2 (23) 



rm (r, „ r,\- (^1-^2)w(P2-P3)mi(P3-Pi)m2 

V,/^2,M3lPl'P2,P3j - —2 ■ (^4j 



Then the three-gluon Green function in the MOM scheme {pi = p\ = p\ = jJ^^) can be 
parametrised as 

AiM)^',M)Kip^)) = r'^[G(^)^^"^(/i^)TW^,,,. (j>i,P2,P3) -Q fs^,^^ _ p_H^\ + 



+ H^'\f^')rif^,,,,M^P2,P3) 



(25) 



The scalar function G^^^^^™{fi^), proportional to the coupling go at tree level, may be 
extracted by the following projection: 



i=l,3 

1 £abc I ^ ^ ^ 



18/i^ 

(26) 
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Asymmetric case 

The three-gluon Green function with one vanishing external propagator (^11,511) can be 
parametrised as 

Gtlf\P. 0, -P) = (aI{-p)AI{pW^{Q)) = 2r'% (^5,. - (27) 
and thus 

1 £abc 

G(3)a.yn.(^2) ^ —, ^^^l^,_^^ S,.P,Gf}f^{p, 0, -p). (28) 

Gauge coupling 

Using the scalar functions G^^^ and G^^^ (the last stands for G'(3)asym Q{3)sym^^ ^j^g gauge 
coupling at the renormalisation scale /i^ is defined by 

n _ G^'^KpIpIpD v3/2. 2^ ..ql 

^«^^^-G(2)(p2)G(2)(p2)G(2)(p|)^3 if^) l^yj 

in the case of three-gluon vertices, where the choice of pi determines the renormalisation 
scheme (MOM or MOM). In paper [213] it is shown that in the MOM schemes gnin'^) 
vanishes in the zero momentum limit. It is interesting to investigate the possible infiuence 
of Gribov copies on this behaviour. We address this question in the following section. 



4.2 Influence of the lattice gauge fixing 

A natural question that arises after the study of the distribution of Fj^in is whether the 
gauge configurations having the same value of Fj^m are equivalent, i.e. they differ only by 
a global gauge transformation. 

This can be checked by calculating the two-point gluonic correlation function on the 
gauge configuration . Indeed, according to the lattice definition of the gauge field that we 
used 

Gflj^ - y) (X Tr [{U,{x) - f/t(x)) ■ (U^y) - Ul{y))] . (30) 

Applying a global gauge transformation U^{x) V^[/^(x)V^ we see that the gluon prop- 
agator remains unchanged. This is also the case of the ghost propagator scalar function. 
We have checked numerically that the values of the gluon and the ghost propagators in 
Fourier space are the same for gauge configurations having the same F-^m (cf. Tables l2f3l 
). For comparison, the data from a lattice Gribov copy (the same gauge configuration but 
having different value of F^i^) is also given. We see that the values are almost the same 
for gauge configurations having the same Fmin contrarily to those coming from different 
minima. In fact, taking in account rounding errors appearing during the calculation, we 
may say that gauge configurations having the same -Fmin are equivalent. 

The second question that has already been considered by ( [^.[8].f9j.[Tfl].[TI].|12|.[T8] 
) is the dependence of Green functions on the choice of the minimum. For this purpose we 
have performed the following simulation: for every of the 100 gauge configurations used 
to compute Green functions the gauge was fixed 100 times, and the Monte-Carlo average 
was computed with respect to the "first copy" (fc) found by the minimisation algorithm 
and the "best copy" (be), having the smallest value of F-^m- We have calculated the gluon 
and the ghost propagators, and also the three-gluon Green functions in symmetric and 
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y •G(2)(p2) : p2 ^ 



p2 = 1 

4 pl41=16 
■ 4 p[41^4 



-0.871010810260 



4249297 

2012518 

1215762 

834876.3 

620065.2 

521585.0 

410698.8 



-0.871010810260 



4249295 

2012516 

1215762 

834876.4 

620067.3 

521585.2 

410698.5 



-0.870645877060 



3322788 

2006186 

1362671 

798032 

434235.6 

556509.6 

440623.9 



Table 2: Gluon correlator calculated on different gauge configurations having the same value of F^\^ 
(columns 2 and 3), compared to values from the "secondary" minimum (column 4). We present data for 
all Hi group orbits for the momentum = A ((2l],C3)- Notice that in all this paper the momenta are 
given in units of 2tt/ (aL). Simulation has been performed with parameters V = W^, (3 ~ 2.4 



F ■ 

^ mm 


-0.871010810260 


-0.871010810260 


-0.870645877060 




p2 = l 


14.06473 


14.06473 


14.82984 




p2 = 2 


6.278253 


6.278253 


6.736338 




p2 = 3 


3.757531 


3.757531 


3.939130 




4 p[41=16 


2.929602 


2.929602 


2.705556 




4 p[4]^4 


2.599088 


2.599088 


2.775566 




p2 =5 


2.071200 


2.071200 


2.011100 



Table 3: Ghost correlator calculated on different gauge configurations having the same value of F^^^^ 
(columns 2 and 3), compared to values from the "secondary" minimum (column 4). We present data for 
all Hi group orbits for the momentum = 4 ((21], 52]). Simulation has been performed with parameters 
V = 16^/3 = 2.4 
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L 






Fi!\p')-Fi'^{p') 


ir(2) 


2.1 


8 




1 


0.211 


0.045 


2.1 


16 


4 


p[41=16 


0.145 


0.033 


2.2 


8 




1 


0.078 


0.019 


2.2 


16 


4 


p[41=16 


0.023 


0.006 


2.3 


8 




1 


0.086 


0.024 


2.3 


16 


4 


pl41 = 16 


0.114 


0.034 



Table 4: Volume dependence of the ghost propagators, from Tables FtI 1261 



asymmetric kinematic configurations. The simulations have been performed on lattices 
of volumes 8^ and 16"^ for /3 = 2.1,2.2,2.3. At these values of j3 we are sure to have 
lattice Gribov copies. The results are given in Tables 1711291 All data is given in lattice 
units. We present data for all H^^ group orbits ([2T],[T3|) for considered momenta, and we 
have checked that our results for two-point Green functions are consistent with numerical 
values given in [Z|,|H|- We conclude from Tables EHSHl that the ghost propagator is quite 
sensitive to the choice of the minimum - in the case of (be) the infrared divergence is 
lessened. No systematic effect could be found for gluonic two- and three-point Green 
functions, the values in the cases of (fc) and (be) being compatible within the statistical 
errors. 

It is conjectured in (22] that in the infinite volume limit the expectation values calcu- 
lated by integration over the Gribov region and the fundamental modular region become 
equal. However, the ghost propagator depends on this choice, even for volumes larger than 
the critical volume defined in section El see Table UHl where there is a four a discrepancy 
for = 1. This dependence has been found to decrease slowly with the volume |T2|. The 
results of section ini indicate that the convergence can only happen beyond the critical size. 
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2.2 8 



-0.8236 
-0.8262 
-0.8272 
-0.8279 



0.003744 
0.002367 
0.001377 
0.000802 



10 
12 
16 



2.4 8 



-0.8642 
-0.8669 
-0.8686 
-0.8702 



0.005270 
0.002739 
0.001849 
0.001003 



12 

12 
16 



Table 5: Volume dependence of the Monte-Carlo+gauge orbit mean value at minima i^min and the 
dispersion of this mean. 

To check this we compare the fc-bc values of the ghost propagator, at one physical value 
of the momentum, for the orbit = 1 on a S'' lattice and the orbit = 4,pW = 16 on the 
16^ lattice ^ at the same /? (see Table IH). It happens indeed that the decrease is observed 
only at P = 2.1 and 2.2, in accordance with Table [H However, these values of P are not 
in the scaling regime, and thus a study on larger lattices would be welcome. It is not 
surprising that the ghost propagator depends on the bc/fc choice: the be corresponds to 
the fields further from the Gribov horizon where the Faddeev-Popov operator has a zero 
mode, whence the inverse Faddeev-Popov operator (ghost propagator) is expected to be 
smaller as observed. The correlation between the bc/fc choice and the gluon propagator 
is not so direct. Another quantity is obviously strongly correlated to the bc/fc choice: the 
value of -Finin. We tested the volume dependence of the Monte-Carlo+gauge orbit mean 
value of the quantity Fmin (see Tab HI). According to the argument given in [22j, all min- 
ima become degenerate in the infinite volume limit, and closer to the absolute minimum 
(in the fundamental modular region). We see from the Table [HI that their average value 
and dispersion decrease with the volume at fixed /?, in agreement with (22] • 

5 Conclusions 

Our study showed that 

• Lattice Gribov copies appear and their number grows very fast when the physical 
size of the lattice exceeds some critical value ~ 2.75/a/o"- This result is fairly 
independent of the lattice spacing. 

• The configurations lying on the same gauge orbit and having the same -Fmin are 
equivalent, up to a global gauge transformation, and yield the same Green functions. 
Those corresponding to minima of Fu with different values of Fmin differ by a non- 
trivial gauge transformation, and thus they are not equivalent. 



• We confirm the result (0, 0,0, [El, dH, [121,111) that the divergence of the ghost 



propagator is lessened when choosing the "best copy" (corresponding to the choice 
of the gauge configuration having the smallest value of Fu). We also showed that 
gluonic Green functions calculated in the "first copy" and "best copy" schemes are 
compatible within the statistical error, no systematic effect was found. We conclude 
that gauge couplings that are defined by amputating the three-gluon vertices only 

^Remember that the momentum in physical units is equal to 27rp/(La) in our notations. 
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slightly depend on the choice of the minimum of Fu- It implies that the infrared 
behaviour of gnifi"^) reported in |2n| is not significantly influenced by lattice Gribov 
copies. 

• We found that the influence of Gribov copies on the ghost propagator decreases with 
the volume when the physical lattice size is larger than the critical length discussed 
above. We also show that the quantity Fmin decreases when the volume increases. 
These two points are in agreement with the argument on the equality of the averages 
over the Gribov's region and the fundamental modular region [22j. 

The fact that the abundance of lattice Gribov copies is mainly an increasing function of 
the physical size of the lattice is not too surprising. The existence of Gribov copies is 
a non-perturbative phenomenon and as such is related |T| to the infrared properties of 
the Faddeev-Popov operator and to the conflnement scale. It should then dominantly 
depend on the infrared cut-off, the size (or volume) of the lattice in physical units. A 
milder dependence on the ultraviolet cut-off, the lattice spacing, is also expected but the 
limited accuracy concerning the lattice spacing did not allow us to identify it. In a recent 
paper the lowest eigenvalues of the Faddeev-Popov operator have been computed. 
A dependence on the size of the lattice is seen. A detailed study of this dependence in 
connection with our flndings in this paper would be useful. 

More generally, the question of extracting from lattice simulations informations about 
Gribov copies in the continuum limit is not a simple issue. However, bulk quantities, 
such as the one we propose in this paper, may be traced down to the continuum limit 
and provide precious information about continuum Gribov copies. The total number 
of lattice Gribov copies in a given gauge orbit may be large. All the minima of the 
functional F_a[u{x)] ^ are possible end-points of the Landau gauge flxing algorithm, but 
their probability to be selected by the algorithm depends on the size of their domain of 
attraction and might vary signiflcantly. This probability is therefore expected to depend 
on the attractor pattern and not on the choice of algorithm, which should be checked. 

The explosion of the number of Gribov copies at larger volume does not contradict the 
statement that they a have decreasing influence on expectation values j22]. Even more, 
this influence decreases when the number of Gribov copies is large enough (cf. Tables lilT|l . 
i.e. above the critical volume that we have outlined. 

We wish to continue this research in a few directions: more reflned quantities may be 
defined, a better accuracy is needed, larger volumes studied, in particular to check if the 
ghost propagatore dependence on the bc/fc choice fades away in the scaling regime. An 
extension of the study to SU{3) and to unquenched gauge configurations will be welcome. 

Note added in proof While this paper was being submitted, a preprint has ap- 
peared [in] which adds to the gauge group the extra global symmetry of the quenched 
lattice action with respect to the centre of the gauge group. This leads to an extension of 
the notion of Gribov copies and it results a sensitivity of the gluon propagator upon the 
latter. We note that the copy dependence is mainly observed at small lattice volumes, 
below the critical size found by us. 
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10- 


-p — 
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10 
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-o 


0.729620 


0.004738 


0.78416564213526 


0.00260692007045 


2.0 


10- 


14 


10- 


5 


10 


-9 


0.729606 


0.004984 


0.78416564421249 


0.00260692062016 


2.0 


10" 


14 


10" 


7 


10 


-9 


0.729606 


0.004985 


0.78416564421259 


0.00260692061867 


2.0 


10- 


18 


10- 


7 


10 


-9 


0.729606 


0.004985 


0.78416564421258 


0.00260692062307 


2.0 


10- 


24 


10- 


7 


10 


-9 


0.729606 


0.004985 


0.78416564421258 


0.00260692062150 


2.0 


10- 


28 


10- 


■7 

I 


10- 


— 9 


0.729606 


0.004985 


0.78416564421258 


0.00260692062134 


2.3 


10- 


10 


10- 


—E 

o 


10 


— p — 

-o 


0.324636 


0.056660 


0.84466892598808 


0.00415444722060 


2.3 


10- 


14 


10- 


5 


10 


-9 


0.324636 


0.056660 


0.84466892602737 


0.00415444721015 


2.3 


10- 


14 


10- 


7 


10 


-9 


0.324636 


0.056660 


0.84466892602833 


0.00415444720619 


2.3 


10- 


18 


10- 


7 


10 


-9 


0.324636 


0.056660 


0.84466892602843 


0.00415444720513 


2.3 


10- 


24 


10- 


7 


10 


-9 


0.324636 


0.056660 


0.84466892602847 


0.00415444720356 


2.3 


10" 


28 


10- 


7 


10 


Q 

— y 


0.324636 


0.056660 


0.84466892602847 


0.00415444720356 


2.8 


10- 


10 


10- 


-p 

O 


10 


— p — 

-o 


0.024946 


0.155566 


0.89393978935534 


0.00777035608940 


2.8 


10- 


14 


10- 


5 


10 


-9 


0.024946 


0.155567 


0.89393978936422 


0.00777035609273 


2.8 


10- 


14 


10- 


7 


10 


-9 


0.024946 


0.155567 


0.89393978936434 


0.00777035609111 


-> 2.8 


10- 


18 


10- 


7 


10 


-9 


0.024946 


0.155567 


0.89393978936448 


0.00777035608963 


2.8 


10- 


24 


10- 


7 


10 


-9 


0.024946 


0.155567 


0.89393978936449 


0.00777035608962 


2.8 


10- 


28 


10- 


7 


10 


-9 


0.024946 


0.155567 


0.89393978936449 


0.00777035608962 



Table 6: The influence of different stopping parameters lfTT|l on the value of the probabiHty P{(5)- 
Simulations done for the lattice size = 8^, 100 Monte-Carlo configurations (we uses the same set of 
configurations for every value of jS) x iV^F = 100 gauge fixings. The last two columns give the average 
value of i^min and the standard dispersion of this average. The arrow indicates our choice of stopping 
parameters ltTT|l . 
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SU{2), y = 8^ and (3 = 2.1, 100 Monte-Carlo configurations x 
gauge fixings 

Two point functions 













1 


4.898 


0.099 


4.687 


0.071 


2 


2.046 


0.039 


1.959 


0.043 


3 


1.210 


0.021 


1.168 


0.023 


4 p[41=16 


0.961 


0.023 


0.925 


0.021 


4 p[41=4 


0.834 


0.019 


0.801 


0.013 


5 


0.696 


0.007 


0.680 


0.014 


Table 7: Ghost propagator, 1/ = 8* /3 = 2.1 















11.161 


0.438 


10.894 


0.418 


1 


6.225 


0.129 


6.248 


0.135 


2 


4.089 


0.035 


4.095 


0.043 


3 


2.883 


0.033 


2.868 


0.023 




2.329 


0.043 


2.305 


0.031 




2.129 


0.023 


2.147 


0.013 


5 


1.773 


0.009 


1.785 


0.008 



Table 8: Gluon propagator, y = 8^ /3 = 2.1 

Three point functions 



p2 










1 


29.636 


2.081 


29.408 


2.721 


2 


19.018 


1.168 


18.387 


1.173 


3 


11.513 


0.971 


12.293 


1.739 


4 p[41=16 


10.424 


0.857 


11.817 


1.126 


4 p[41^4 


6.968 


0.618 


6.547101 


0.471 


5 


5.164 


0.264 


4.697285 


0.293 



Table 9: Three-gluon vertex in asymmetric kinematic configuration, y = 8"* /? = 2.1 





G(3)Bym^p2) 








2 


13.613 


0.745 


14.188 


0.752 


4 


2.680 


0.272 


2.612 


0.272 



Table 10: Three-gluon vertex in symmetric kinematic configuration, F = 8* /3 = 2.1 
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SU{2), y = 8^ and (3 = 2.2, 100 Monte-Carlo configurations x 100 
gauge fixings 

Two point functions 





4:\p') 






Ki\p') 


1 


4.178 


0.103 


4.100 


0.070 


2 


1.749 


0.045 


1.735 


0.037 


3 


1.042 


0.027 


1.032 


0.024 


4 p[41=16 


0.835 


0.030 


0.821 


0.024 


4 p[41=4 


0.720 


0.016 


0.715 


0.014 


5 


0.614 


0.017 


0.609 


0.014 


Table 11: Ghost propagator, V = 8^ f3 = 2.2 




G^fip') 


SG^fip') 


gZHp') 


SG'^^P') 





20.170 


0.620 


20.241 


0.665 


1 


8.359 


0.129 


8.247 


0.122 


2 


4.431 


0.039 


4.433 


0.037 


3 


2.870 


0.028 


2.857 


0.025 




2.111 


0.038 


2.153 


0.041 




2.021 


0.013 


2.038 


0.019 


5 


1.642 


0.014 


1.652 


0.013 



Table 12: Gluon propagator, V = 8^ P = 2.2 
Three point functions 



p' 










1 


78.913 


4.405 


81.523 


5.131 


2 


36.170 


1.393 


35.784 


1.809 


3 


17.851 


2.147 


18.475 


1.860 


4 p[41=16 


20.516 


1.064 


20.621 


1.115 


4 p[41^4 


13.097 


0.645 


13.382 


0.595 


5 


8.920 


0.336 


8.985 


0.383 



Table 13: Three-gluon vertex in asymmetric kinematic configuration, V — 8'^ P — 2.2 



p' 


G(3)Bym^p2) 








2 


21.099 


1.212 


20.483 


1.148 


4 


2.494 


0.185 


2.512 


0.189 



Table 14: Three-gluon vertex in symmetric kinematic configuration, V = 8'^ f3 = 2.2 
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SU{2), y = 8^ and (3 = 2.3, 100 Monte-Carlo configurations x 100 
gauge fixings 

Two point functions 











Ki\p') 


1 


3.742 


0.109 


3.656 


0.084 


2 


1.569 


0.043 


1.544 


0.032 


3 


0.945 


0.022 


0.941 


0.019 


4 p[41=16 


0.755 


0.021 


0.751 


0.016 


4 p[41=4 


0.669 


0.015 


0.663 


0.013 


5 


0.562 


0.009 


0.564 


0.009 


Table 15: Ghost propagator, y = 8^ /3 = 2.3 















35.140 


1.555 


34.659 


1.597 


1 


9.7073 


0.107 


9.839 


0.156 


2 


4.487 


0.054 


4.503 


0.055 


3 


2.668 


0.027 


2.652 


0.018 




1.934 


0.022 


1.919 


0.020 




1.811 


O.OK) 


1.823 


0.013 


5 


1.414 


0.014 


1.414 


0.014 



Table 16: Gluon propagator, V = 8^ P = 2.3 
Three point functions 



9 

P 










1 


614.254 


42.333 


641.328 


47.008 


2 


165.962 


7.456 


161.354 


7.730 


3 


65.054 


2.786 


63.734 


2.734 


4 p[41=16 


33.456 


3.787 


32.685 


3.535 


4 p[41^4 


33.770 


1.801 


35.883 


2.099 


5 


20.471 


0.9744 


19.994 


1.168 



Table 17: Three-gluon vertex in asymmetric kinematic configuration, V — 8"^ (3 — 2.3 



p' 


G(3)Bym^p2) 








2 


26.483 


1.260 


26.774 


1.280 


4 


2.249 


0.105 


2.424 


0.104 



Table 18: Three-gluon vertex in symmetric kinematic configuration, V = 8'^ f3 = 2.3 
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SU{2), V = 16^ and (3 = 2.1, 100 Monte-Carlo configurations x 
gauge fixings 

Two point functions 











Ki\p') 


1 


23.930 


0.304 


22.615 


0.226 


2 


10.538 


0.200 


10.188 


0.113 


3 


6.372 


0.159 


6.257 


0.087 


4 p[41=16 


4.551 


0.116 


4.405 


0.103 


4 p[41=4 


4.489 


0.119 


4.421 


0.062 


5 


3.437 


0.086 


3.344 


0.053 


Table 19: Ghost propagator, y = 16^ /3 = 2.1 


9 













8.485 


0.440 


8.309 


0.251 


1 


8.182 


0.122 


7.871 


0.102 


2 


7.186 


0.073 


7.022 


0.071 


3 


6.359 


0.037 


6.385 


0.037 


4 11=10 


6.058 


0.091 


5.874 


0.086 




5.710 


0.081 


5.G()9 


0.007 


5 


5.161 


0.029 


5.195 


0.026 



Table 20: Gluon propagator, y = 16^ /? = 2.1 
Three point functions 



p' 










2 


19.861 


3.809 


27.689 


5.911 


4 


14.769 


2.586 


18.952 


2.937 



Table 21: Three-gluon vertex in symmetric kinematic configuration, V = 16* /3 = 2.1 
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SU{2), V = 16^ and (3 = 2.2, 100 Monte-Carlo configurations x 
gauge fixings 

Two point functions 











Ki\p') 


1 


20.867 


0.376 


20.406 


0.392 


2 


9.151 


0.145 


9.070 


0.192 


3 


5.499 


0.083 


5.460 


0.085 


4 p[41=16 


3.861 


0.095 


3.838 


0.082 


4 p[41=4 


3.803 


0.067 


3.819 


0.068 


5 


2.929 


0.041 


2.979 


0.067 


Table 22: Ghost propagator, y = 16^ /3 = 2.2 






SG^fip') 









14.473 


0.676 


15.380 


0.635 


1 


12.614 


0.255 


12.330 


0.184 


2 


10.564 


0.066 


10.531 


0.097 


3 


8.813 


0.081 


8.769 


0.061 




7.760 


0.162 


7.577 


0.089 




7.447 


0.032 


7.429 


0.073 


5 


6.393 


0.052 


6.395 


0.044 



Table 23: Gluon propagator, y = 16^ /? = 2.2 
Three point functions 



p' 










1 


136.266 


91.211 


115.376 


53.064 


2 


93.912 


39.689 


90.432 


26.263 


3 


84.119 


17.125 


81.880 


17.916 


4 p[41=16 


101.745 


17.120 


99.140 


18.988 


4 p[41^4 


78.176 


15.640 


63.382 


14.634 


5 


51.954 


7.666 


56.088 


8.017 



Table 24: Three-gluon vertex in asymmetric kinematic configuration, V = 16^ /3 = 2.2 



p' 


G(3)Bym^p2) 








2 


93.161 


11.297 


87.022 


10.680 


4 


34.756 


7.434 


37.655 


5.797 



Table 25: Three-gluon vertex in symmetric kinematic configuration, V = 16^ /3 = 2.2 
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SU{2), V = 16^ and (3 = 2.3, 100 Monte-Carlo configurations x 
gauge fixings 

Two point functions 











Ki\p') 


1 


18.326 


0.297 


17.157 


0.172 


2 


7.993 


0.150 


7.645 


0.089 


3 


4.824 


0.086 


4.684 


0.060 


4 p[41=16 


3.431 


0.053 


3.317 


0.052 


4 p[41=4 


3.364 


0.054 


3.338 


0.048 


5 


2.581 


0.025 


2.533 


0.036 


Table 26: Ghost propagator, y = 16^ /3 = 2.3 






SG^fip') 









32.705 


1.063 


31.557 


0.721 


1 


21.940 


0.444 


21.616 


0.385 


2 


15.007 


0.122 


15.206 


0.166 


3 


11.400 


0.107 


11.235 


0.116 




8.848 


0.190 


8.882 


0.108 




8.904 


0.118 


8.782 


0.099 


5 


7.057 


0.062 


6.978 


0.045 



Table 27: Gluon propagator, y = 16^ /? = 2.3 
Three point functions 



p' 










1 


1192.828 


250.675 


779.501 


204.768 


2 


566.632 


39.161 


505.435 


46.025 


3 


402.454 


32.415 


388.260 


26.754 


4 p[41=16 


272.205 


35.491 


242.150 


30.928 


4 p[41^4 


239.204 


35.152 


241.295 


25.467 


5 


174.210 


14.126 


169.739 


10.825 


6 


130.209 


7.523 


125.733 


7.563 



Table 28: Three-gluon vertex in asymmetric kinematic configuration, V = 16* (3 = 2.3 



p' 










2 


283.901 


15.417 


313.824 


15.595 


4 


81.695 


5.785 


72.573 


5.787 



Table 29: Three-gluon vertex in symmetric kinematic configuration, V = 16* (3 = 2.3 
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